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I.  INTRODUCTION 


During  Che  spring  o£  1983,  the  Ballistic  Research  Laboratory  measured 
behind-the-armor  effects  for  several  weapons  against  a  variety  of  targets. 

One  of  the  major  concerns  was  heat  generation  during  and  immediately  following 
impact  of  the  warhead  against  the  targets. 

Previous  tests  of  this  type  had  used  thin  skin  total  heat  gage9  and 
standard  laboratory  heat  detectors.  Thin  skin  gages  record  only  the  maximum 
temperature  attained  by  skins'  back  surface  with  no  indication  of  the  time 
reauired  to  reach  this  temperature.  Although  an  approximation  of  the  total 
heac  deposition  can  be  made,  no  rate  of  deposition  can  be  estimated  without 
knowing  the  time  of  the  event.  The  large  thermal  mas9  of  a  typical  laboratory 
detector  such  as  an  infrared  power  meter  is  designed  to  measure  a  steady-state 
flux  and  cannot  respond  quickly  enough  to  register  a  transient  event.  Neither 
gage  type  can  withstand  severe  blast  and  shock  environments. 

'  i  '  ■ 

Since  the  temperature  inside  the  target  rose  sharply  in  a  very  9hort 
time,  special. gages  were  fabricated  which  emphasized  the  time  response  of  the 
gage  rather  than  its  thermal  capacity.  This  gage  is  used  when  the  thermal 
flux  is  of  short  duration,  150  milliseconds  or  les9,  and  of  low  total  energy, 
179  Calories  per  square  centimeter  or  less.  The  maximum  operating  temperature 

of  the  gage  is  500  degrees  Celsius.  ' 

"  ''  •  •  *  r,.  ^  .  » 

II.  GAGE  DESCRIPTION 


A.  Basic  Gage  Configuration 


The  basic  configuration  of  the  thermal  conductivity  gage  is  shown  in 
Figure  1.  The  measuring  element  of  the  gage  is  a  thin  aluminum  disc  with  a 


Figure  1.  Basic  Thermal  Conduct ivity ' Cage  Configuration 
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chromel-alumel  thermocouple  mounted  across  the  center  of  each  face  using  a 
thermoconducting  epoxy.  This  subassembly  is  mounted  on  the  end  of  a  ceramic 
mandrel  which  is  inserted  into  the  ceramic  sleeve  so  that  the  front  surface  of 
the  disc  is  ti>.3h  with  the  front  of  the  sleeve.  Wires  connecting  the  gage  to 
any  external  ir strumentation  are  anchored  to  the  mandrel  and  then  passed  out 
Che  rear  of  the  sleeve.  Detailed  drawings  of  the  gage  are  in  the  Appendix. 

B.  Design  Calculations 

Eased  on  the  following  assumptions  and  material  properties,  gage 
performance  was  calculated  from  the  general  one-dimensional  thermal 
conductivity  equation, 

3T/3 1  -  a32T/3x2  (1) 

where 

T  *  temperature  in  degrees  Celsius; 
t  *  time  in  seconds; 

a  •  thermal  diffusivity  in  square  centimeters  per  second;  and 
x  *  distance  in  centimeters. 


The  equation  was  solved  by  the  graphical  method  described  in  Reference  1  and 
illustrated  in  Figure  2.  Time  response  of  the  gage  was  calculated  by  plotting 
temperature  distribution  versus  time  as  heat  passes  longitudinally  through  the 
different  materials  in  the  gage  in  constant  differential  increments  of  time. 
This  differential  time  element  and  the  incremental  distance  through  the  gage 
are  related  by 


X  m  mm 

Figure  2.  Thermal  Geometry  for  Calculation 


M.  Jakob,  G.E.  Hawkins,  Elements  of  Heat  Transfer,  John  Wiley 
&  Sons,  Inc.,  New  York,  1958. 


1.  Assumptions.  The  following  assumptions  were  u9ed  in  the  calculations. 

a.  The  heat  flow  in  and  out  of  the  gage  ia  parallel  to  the 
longitudinal  axis  of  the  gage.  Heat  flow  in  the  radial  direction  is 
ins ignif icant . 


b.  The  cement  used  to  bond  the  thermocouples  to  the  disc  has 
conductive  properties  similar  to  aluminum.  The  thickness  of  the  bond  is 
approximately  one  wire  diameter. 

c.  The  temperature  measured  by  the  thermocouples  is  assumed  to  be 
the  average  temperature  in  the  bond  layers. 

d.  The  primary  sources  of  heat  are  radiant  in  nature,  so  film  and 
boundary  layer  effects  are  ignored  in  calculating  the  thermal  response 
characteristics.  These  effects  retard  convective  and  conductive  heat  transfer 
but  do  not  significantly  influence  radiant  heat  transfer. 

e.  The  thermocouple  wires  are  of  sufficiently  small  cross  section 
and  of  sufficient  length  to  render  heat  loss  from  the  gage  insignificant 
(0.076  millimeters  in  diameter  ana  30  millimeters  long). 

2.  Material  Properties.  The  disc  is  manufactured  from  2024-T4  aluminum 
alloy.  Both  the  mounting  mandrel  and  the  surrounding  sleeve  were  fabricated 
from  Macor,  an  alumina-baaed,  machinable  ceramic.  Denaity,  thermal 
conductivity,  and  specific  heat  were  obtained  from  the  cited  literature. 
Thermal  diffusivity  is  determined  by  the  relation 

a  -  k/pC  .  (3) 

P 

The  2024-T4  aluminum  properties  are:^ 


density  (p) 
thermal  conductivity  (k) 
specific  heat  (Cp) 
thermal  diffusivity  (a) 

•j 

The  Macor  properties  are. 


0.135  gm/cra  ; 

2.64  Cal/cm  sec*C; 
32.8  Cal/gm*C;  and 
0.597  cm  /sec. 


2 

P  ■  0.125  gm/cm  ; 
k  ■  3.1  x  10"3  Cal/cm  aec'C; 

C  -  180  to  200  Cal/gm‘0;  and 
a  ■  0.11  to  0.12  x  10-3  cm^/sec. 


t) .  A .  Maykuth,  Aerospace  Structural  Metals  Handbook,  Bel  four 
atuien,  Inc.,  Columbus,  OH,  1  9  8  0  , 

Corning  Bulletin,  "Macor  Machinable  Glass  Ceramic,"  Corning 
Glass  Works,  Corning,  NY,  1981. 
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The  disc  is  0.0305  centimeters  wide  and  has  a  diameter  of  0.625 
centimeters.  The  thermocouple  wire  diameter  and  the  bond  layer  thickness  are 
both  0,0076  centimeters. 

3.  Response  Calculation  Results.  The  geometry  in  Figure  2  and  the  above 
stated  material  properties  were  used  in  the  graphical  solution  of  Eq.  (1)  for 
the  conditions  of  changing  surface  temperature  and  conduction  layers  of 
ciffering  materials.  Figure  3  ehows  the  response  of  the  front  thermocouple  to 
a  surface  temperature  step  rise  of  one  millisecond  and  Figure  4  shows  the 
response  of  the  same  thermocouple  to  a  250  Herts  sinusoidal  surface 
temperature  variation.  The  shaded  portiona  of  the  curves  show  the  uncertainty 
in  the  calculations  due  to  variations  of  the  actual  material  properties  and 
vagaries  introduced  by  the  assumptions. 

4.  Calculation  of  Flux  Penalty.  The  flux  density  (q')  is  the  sura  of  the 
instantaneous  conductive  heat  flow  density  (q‘)  and  the  differential  heat 
storage  density  (q^)  •  These  quantities  are  defined  by  the  following 
equations: 

2 

q'  ■  qj  +  Cal/cm  sec;  (4) 

qjj  -  kAT/Ax  Cal/cra2aec;  (5) 

q'  -  (pc  /2)(AT  /At  +  AT,/At)Ax  Cal/cm2sec;  (6) 

where  T^  is  the  front  thermocouple  temperature  and  T^  is  the  rear  thermocouple 
temperature.  The  linearity  shown  in  the  above  equations  is  true  for  steady- 
state  heat  flow  and  for  sufficiently  small  values  of  Ax.  By  combining  Eqs. 
(4),  (5),  and  (6),  the  following  equation  is  derived  for  total  flux  density: 

q'  -  ( k/Ax) (T.  -  T.)  ♦  (pC  Ax/2)(AT,/At  ♦  AT  /At)  Cal/cm2sec;  (7) 

lip  i  2 

where  Ax  (■  0.044  cm)  is  the  distance  between  the  two  thermocouples . 
Evaluating  the  constants  in  Eq.  (7)  yields 

q'  -  SO.OCTj  -  T2)  ♦  0.097(AT1/At  f  AT2/Ac)  Cal/cm2sec.  (8) 

5.  Calculation  of  Total  Deposition  Pen8ity.  The  total  heat  deposition  (O')  is 
the  time  integral  of  the  flux  density  ( q  * } . 
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where  (t^  -  tQ)  is  the  time  interval  of  observation.  This  integration  yields 


i 

r. 

i 


Q'  -  60.0/tf(T1  -  T2)dt  ♦  0.097(T1  +  Tj) 


*  Cal/cm2. 


(10) 


The  interval  of  interest  for  both  thermocouples  is  the  time  between  the  first 
rise  from  the  baseline  temperature  (tQ)  and  the  point  at  which  the  maximum 
temperature  is  reached  (t£)  by  the  front  thermocouple.  Since  the  integration 
rero  time  is  taken  at  the  start  of  the  event,  reference  junctions  in  the  data 
acquisition  system  are  unnecessary  as  long  as  the  changes  in  temperature  are 
referenced  with  respect  to  the  temperature  at  time  tQ. 

6.  Approximate  Solutions.  Since  the  measurand  is  a  transient  and  the  aft 
portion  of  the  gage  has  a  much  lower  thermal  conductivity  than  the  front 
portion  of  the  gage,  the  deposition  and  flux  density  are  actually  nonlinear 
and  of  the  forms: 


q'  -  f{(Tj  -  I2),  (AT^/At  ♦  AT^/At) }  Cal/cm2sec, 

and 


(11) 


k 


O'  -G{  /  f(T:  -  T2)dt,  (Tj+  T2>}  Cal/ cm‘ 


(12) 


The  existence  of  significant  levels  of  interference  in  the  data  acquired  from 
the  field  make  Eqs.  (11)  and  (12)  quice  cumbersome  to  evaluate.  Simpler 
algorithms  can  be  developed  based  on  the  limits  of  Eqs.  (11)  and  (12)  in  the 
time  domain  and  least  squares  fits  to  calibration  data. 

From  Eq.  (11),  it  becomes  evident  that  for  short  time  intervals,  the 
average  temperature  in  the  gage  becomer  the  dominant  parameter.  The  flux 
density  can  be  expressed  as: 
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sec 


(13) 


q'  -  f{  (ATj/At  +  AT2/At)/2j  Cal/cm2 
and  Che  total  deposicion  as: 

Q'  -  G{(ATi  •*-  &T2)/2}  Cal/cm2,  (14) 

where 

AT  *  T  -  T  over  the  time  interval  At; 

*  *  a  + 

f  o 

AT2  ■  T2  ~  t2  ovcr  c^e  time  interval  At;  and 
r  o 

At  ■  t-  -  t  ,  the  time  interval, 
f  o 

These  approximations  resemble  the  equation j  for  thin  skin  gages,  except 
for  the  introduction  here  of  a  time  variable. 


III.  GAGE  CALIBRATION 


A.  Calibration  Procedure 

• 

The  gage  was  calibrated  by  comparing  its  output  to  the  radiant  energy 
output  of  a  carbon  dioxide  laser  operating  at  10.6  micrometers.  A  fast- 
response  photon  drag  detector  was  used  to  measure  the  laser  energy  directly; 
the  laser  energy  was  also  calculated  from  the  laser  output  power  and  th?  time 
duration  of  the  pulse.  The  power  is  constant  for  any  one  calibration  frial 
and  the  pulse  time  is  equivalent  to  the  current  input  time.  The  gage  outputs 
were  recorded  through  data  amplifiers  into  a  transient  recording 
oscilloscope.  The  recording  components  were  calibrated  against  a  digital 
voltage  standard  both  individually  and  as  a  system.  Figure  S  shows  the 
apparatus  used  in  the  calibration. 

B*  Apparatus 

The  following  is  the  list  of  equipment  used  in  the  calibration  procedure. 

1.  CC>2  Laser:  Apollo  Lasers  Inc.,  Model  300,  serial  number  101. 

2.  Infrared  Power  Meter:  Coherent  Rsdiacion  Laborului ics  Inc.,  detector 
serial  number  081. 
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3.  Germanium  Focusing  Lens:  Coherent  Radiation  Laboratories,  Inc., 

10  inch  focal  length,  no  serial  number. 

4.  Transient  Recording  Oscilloscope:  Nicolet  Corp.,  Model  Explorer  II, 
serial  number  612831. 

5.  Instrumentation  Amplifiers:  Dynamics  Inc.  Model  7525,  serial  numbers 
3543  and  3563. 

6.  DC  Voltage  Standard:  Analogic  Inc.,  Model  AN3100,  serial  number 
7825476. 

7.  Photon  Drag  Detector:  Rotin,  Inc.,  Model  741D,  aerial  number  081. 


CCRUANIUU 


C.  Recording  System  Calibration 

The  recording  system  was  calibrated  on  site,  in  three  separate  steps. 

First,  the  recording  oscilloscope  was  calibrated  against  the  voltage 
standard  on  both  channels  A  and  B  over  all  ranges.  These  calibrations, 
performed  in  accordance  with  the  manufacturer's  instructions,  included  both 
the  recording  and  the  playback  systems.  This  procedure  must  be  completed 
before  calibrating  the  instrumentation  amplifiers. 

Next,  the  instrumentation  amplifiers  were  calibrated  for  zero,  balance, 
linearity,  and  gain  using  the  manufacturer's  procedures  as  outlined  in  the 
equipment  manual.  The  voltage  standard  was  used  as  the  reference. 
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Finally,  the  acquisition  and  recording  system  was  assembled  as  it  was  to 
be  used  during  the  testa  and  checked  for  ground  loop  and  noise  interference. 

A  ten-step  voltage  pulse  was  recorded  on  every  track  of  each  channel.  The 
pulses  were  stored  on  floppy  disc  and  later  used  to  establish  calibration  and 
conversion  coefficients  for  data  analysis. 

Figure  6  shows  the  recording  system  calibration  scheme. 


Figure  6.  Recording  System  Calibration  Scheme 


The  overall  system  accuracy  was: 

Linearity  -  0.2Z  full  scale; 

Amplitude  Resolution  -  0.05X  full  scale; 
Time  Resolution  -  0.2Z  range; 

Amplifier  Gain  -  100  +  0.02X. 


D.  Laser  Calibration 


Since  the  laser  output  cannot  be  monitored  without  disturbing  the  beam, 
the  laser  was  calibrated  by  determining  the  relationship  between  output 
radiation  and  input  current.  The  photon  drag  detector  provided  the  output 
data  and  the  laser  control  panel's  current  monitor  supplied  both  the  input 
data  and  the  triggering  pulse  to  the  recording  system.  In  both  the  single 
pulse  and  the  multipulse  modes,  the  duration  of  the  output  pulse  was  about 
250  microseconds  less  chan  that  of  the  current  pulse.  The  output  pulse  showed 
a  delay  of  about  50  microseconds  from  the  input  pulse. 

For  each  calibration  test  series,  the  laser  beam  was  focused  over  the 
active  area  of  the  temperature  gage.  The  gage  was  then  removed  and  the  beam 
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placed  in  multipulse  mode.  The  energy  delivered  by  each  pulse  was  calculated 
using  the  average  power,  the  frequency,  and  the  period  of  the  duty  cycle.  The 
power  was  measured  by  the  infrared  meter  and  the  last  two  quantities 
determined  by  the  photon  drag  detector. 

Figure  7  shows  a  schematic  of  the  calibration  setup. 


IR  POWER  |R  PHOTO- 


Figure  7.  Laser  System  Calibration  Scheme 


B.  Cage  Calibration 

Once  the  laser  output  parameters  were  established,  tlje  test  gage  was 
inserted  into  the  beam  as  shown  in  Figure  8.  The  gage  output  voltages  were 
recorded  for  both  single  and  multipulse  modes.  Figures  9  and  10  show  the 
typical  gage  response  to  a  single  pulse  and  to  multiple  pulses,  respectively. 

The  results  of  the  calibration  trials  are  shown  in  Figure  11. 
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IV.  RESULTS 


A.  Gage  Response 

The  operating  response  of  the  gage  ranged  from  one  to  250  milliseconds. 
Pulses  shorter  than  one  millisecond  could  not  be  measured,  and  pulses  longer 
than  250  milliseconds  caused  saturation  of  the  gage,  with  full  saturation 
occurring  at  about  200  milliseconds. 

Figures  12  and  13  show  typical  field  measurements  for  the  front  and  rear 
thermocouples.  The  calculation  shown  in  Figures  3  and  4  indicates  that  the 
actual  peak  surface  temperature  of  the  gage  can  be  20*  higher  than  the  front 
thermocouple  measures.  This  discrepancy  is  larger  during  the  pulse's  rise. 
Because  of  this  uncertainty  in  the  measured  temperatures,  the  approximations 
for  hest  flux  density  and  heat  deposition  yielded  only  average  values  over  the 
interval  At  shown  in  Figures  12  and  13. 

B.  Evaluation  of  Peat  Flux  Density 

The  approximation  for  flux  density,  Eq.  (13),  and  the  analytical 
relation,  Eq.  (7),  do  not  describe  the  gage  operation  well.  This  is  due  to 
the  extreme  sensitivity  of  differentiation  to  signal  noise  and  to 
discontinuities  caused  by  digitization.  A  satisfactory  estimate  of  the 
average  flux  density  was  obtained  by  dividing  the  total  deposition  by  the  time 
of  deposition,  At: 


q'  ■  Q' / At  Cal/cm  sec 


(15) 


Figure  12.  Typical  Field  Measurement  by  Front  Thermocouple 
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C.  Evaluation  of  Heat  Deposition 


The  approximation  for  heat  deposition,  Eq.  (14),  provided  a  reasonable 
description  of  the  gage  operation*  Least  squares  calculations  showed  the 
deposition  density  function  and  its  coefficients  to  be 


<,• .  0.0,18  -.y»s*  cue2, 


(16) 


where  ATj  and  AT2  were  determined  as  shown  in  Figures  12  and  13.  Figure  14 
compares  this  calculated  deposition  density  with  the  measured  data. 


V.  CONCLUSIONS 

The  gage  as  presently  constituted  performs  reasonably  well  as  a  thermal 
flux  gage,  although  design  changes  could  simplify  both  the  data  acquisition 
and  analysis.  Such  changes  include  adoption  of  a  two-wire  shielded 
thermocouple  and  vapor-deposited  thermocouple  junctions.  The  two-wire  system 
would  reduce  the  number  of  necessary  data  channels  to  one  per  gage,  thus 
eliminating  the  differencing  step  in  the  analysis.  Use  of  the  vapor-deposited 
junctions  would  allow  more  accurate  measurement  of  the  gage  surface 
temperature. 
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APPENDIX 

Figures  A1  through  A5  are  detailed  engineering  drawings  of  the  gage  fully 
assembled  and  its  constituent  parts. 
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Figure  A-3.  Thermocouple  Hiring  Mandrel 
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